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Abstract 

An analytically expression for the temperature dependence of the signal gain of an 
erbium-doped fiber amplifier (EDFA) pumped at 1480 nm are theoretically obtained by 
solving the propagation equations with the amplified spontaneous emission (ASE). It is 
seen that the temperature dependence of the gain strongly depends on the distribution 
of population of Er 3+ -ions in the second level. In addition, the output pump power 
and the intrinsic saturation power of the signal beam are obtained as a function of the 
temperature. Numerical calculations are carried out for the temperature range from 
— 20 to + 60 °C and the various fiber lengths. But the other gain parameters, such as 
the pump and signal wavelengths and their powers, are taken as constants. It is shown 
that the gain decreases with increasing temperature within the range of L < 27 m. 
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1 Introduction 

EDFAs have been presenting many advantages such as high gain and low noise in the op- 
tical communication networks and providing a broadband amplification of radiation whose 

wavelength is in the so-called third window for fiber-optic communication (~1530 nm). In 
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addition, the temperature dependence of the gain characteristics of EDFAs has also of great 
importance for WDM systems |IJ. An analytical solution of the rate equation has been 
derived and the gains at optimum amplifier lengths have been experimentally found for the 
various temperature values in the previous works [21 IH E] . The theoretical and experimen- 
tal results for the temperature dependence of the gain with the various lengths of EDFAs 
have been reported |§J, but the temperature dependent analytic expressions has not been 
given in terms of the Boltzmann factor. Afterwards, the theoretical analysis of amplification 
characteristics of EDFAs has been developed to calculate the signal gain, using the rate 
equation model, [7|, and this rate equation has been modified by including the temperature 
and cross section factors to understand the dependence of the gain on the temperature (HI Ej 
for EDFAs pumped at 1480 nm. In this article, we present an analytical expression for the 
signal gain in EDFAs, using the propagation equations improved by including the temper- 
ature effects, and the numerical results for the temperature ranges of — 20 °C to +60 °C. 
We took into account the amplified spontaneous emission (ASE), but neglected the excited 
state absorption (ESA) effect for the simplicity 



The simplest treatment of EDFA considers the two-level amplification system with energy 
levels as shown in Figure 1, when it is pumped at 1480 nm. In this figure, level 1 is the 
ground level and level 2 is the metastable level characterized by a long lifetime r (= 7 _1 ), 
Rp ,e is the pump absorption and stimulated emission rates, 612,21 is the signal stimulated 
absorption and emission rates, respectively. N 2 + and N 2 - are the populations of Er 3+ ions 
within the sub-levels of the second energy state and it is possible to consider each of them as 
a single energy level. Actually, this system contains many sub-levels where the erbium-ions 
reside and they are unequally populated due to the thermal distribution of the ions. Thus, 
the relative occupation of the sub-levels in the thermal equilibrium must be arranged as a 
function of the temperature. This arrangement is governed by Boltzmann's distribution law: 



2 Theory 
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where T is the temperature in degrees Kelvin, fee is Boltzmann's constant. E 2+ and 
E2- are the higher and lower sub-levels energies of the second level, respectively, and 
AE 2 = E 2+ — E 2 _ [10J. C+ r and C~ r are the nonradiative rates which correspond to the 
thermalization process occurring within each manifold of the second level. The rate equa- 
tions corresponding to the two levels 1 and 2 can be given as follow 

^± = r;n x - r;n 2+ + c- iv 2 _ - c+n 2+ , (2) 

" ' A ' S 12 N X - S 21 N 2 _ - iV 2 _ 7 - C~ N 2 _ + C+ iV 2+ , (3) 



dt p ^ p 
Thus, at stationary conditions we now obtain 



R;N 2+ - R"Nx + S 21 N 2 _ - S 12 N! + iV 2 _ 7 . (4) 
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where the populations are time invariant, i.e., dNi/dt — (i — 1,2). In the last two equa- 
tions, bp ,e = hi/p/T<jp ,e , bf e = hu s /rag' e , v p and v s are the pump and the signal frequencies, 
respectively; a^ e is the stimulated absorption and emission cross sections of the pump beam 
while <r" ,e is the the stimulated absorption and emission cross sections of the signal beam, 
respectively; I p and I s are the pump and signal intensities and I ase * s the forward (+ sign) 
and backward (- sign) propagating optical intensities, respectively, rj is the ratio between 
the signal emission and absorption cross sections, and the total concentration distribution 
of Er 3+ ions is N, N = N x + iV 2 „ + N 2+ or in terms of /3, N = JVi + (1 + /3)JV 2 _. 

The differential equations for propagation of the signal, pump and ASE powers are given, 
respectively, as follows 

d P r°° 

= 2tt j o I. [oZN*-. (r) - a a s Nt (r )] rdr, (7) 

d P r°° r 

= ±2tt J o I p [(3a e p N 2 4r) - a a p N x {r) 
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rdr, 



dp 



ASE 



dz 



/■OO poo 

±2hu s 2ira e s N 2 _fi SE (r)rdr±2ir / KiV 2 _(r) - o^r)] P^r*, (9) 
jo Jo 



where f^sE * s t ne normalized ASE intensity profile, P^se ls the amplified spontaneous 
emission power at the position z and has to be determined from a forward as well as a 
backward travelling ASE spectrum, 



P ± 

ASE 



Pase + P. 
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We can decompose the intensity as I SyP (z,r) = P SjP (z)f S}P (r) where P s ^ p (z) is z-dependent 
signal od pump powers and f StP {r) is the normalized signal and pump transverse intensity 
profiles, respectively At this point, by substituting N — N% + (1 + 0)N 2 - into Eq.((7J), we 
have the propagation equation for the signal power: 

dP a 



dz 



poo 

2na a s P s (l + {3 + ri) N 2 „f s (r)rdr - P s 
Jo 



0:.s 



(11) 



where a s = 2-xo a s Jq 00 N(r)f(r)rdr is the absorption constant of the signal beam. To evaluate 
the integral at the right-hand side of Eq. (jllj) . we make use of Eq.(j3J). In this case, multiplying 
both-hand side of Eq.fjSJ) with rdr and then integrating between and oo, we obtain the 
following equations: 



poo poo q- 1 pc 

/ N 2 _rdr = / -r-Z-KN! - (3a e p N 2 _)rdr + / t^K^i - <N 2 _)rdr 
Jo Jo hv v F F Jo ti^s 



+ 



oo T J 
oo T J 



ASE i a 



o hv s 



<iVx - a e s N 2 _)rdr, 



(12) 
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T dP p T dP s T dP^SE 
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/ N 2 ^f(r)rdr = , 

lo Jy ' 2it(A/T -2ra\ 



1 dP p + J_fdP 1+ dPZ SE \ 



hv p dz 



hv K \ dz 



dz J 



(14) 



where we define the confinement factor T = A J °° N 2 _f(r)rdr/ / °° N 2 _rdr and A is the 
effective doped area. We can put the equations into more practical form supposing the 
pump, signal and ASE profiles to be approximately equal, so that the transverse profiles 
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fpi r ) ~ fs{ r ) ~ fASE( r ) = f( r ) an d considering the co-propagating scheme in the positive z 
direction for the simplicity. Inserting Eq. (jl4)) into Eq. fjlljl . we have 



dPs 
dz 



-P, U s + 



hv s 

pint 



1 dP p 
hv p dz 



1 fdP s + dP+ SE 



hv s \ dz 



dz 



(15) 



where the intrinsic saturation power of the signal beam is introduced as follow 



yillt 



hv s {A - 2r(T s T)/r<r(l + (3 + 77). 



(16) 



Therefore, we define the intrinsic saturation power as a function of the temperature. Integrat- 
ing Eq. (J15|) . we obtain the output signal power at z = L and hence establish a relationship 
between the amplifier gain and length: 



int 



P p (0) - P P (L) + (P s (0)+P+ 5g (0)) - (P S (L) + PX SE {L)) 



hv„ 



(17) 



The amplifier gain G = P s (L)/P s (0) can be calculated from the following equation: 



G = exp(—a s L)exp 



with boundary condition 



hu. 



P. 



int 



Pp(0)-P p (L) P 8 (Q) 



l-G) 



Pase(L) 



ASE 



(0) = 0. 



(19) 



If one neglects the effect of the f3 parameter and the ASE power in the gain equation, it can be 
easily seen that the relevant equation is reduced to the previous works Thus, Eq.(fT9"]) 

is a more accurate solution for the propagation equations. In order to obtain the output 
pump power P p {L) in Eq. ljTSj) for the maximal pumping efficiency, it should be substituted 
Eq.® into Eq.flTJ) and Eq.©, and then Eq.© divided by Eq.©. If the obtained result 
makes equal to zero, we have 

P P (L) = , 1 — , (20) 



R ( -a- - £ 



where R = N(r)f(r)rdr/ J °° N(r)rdr. It is notes that the output pump power is a 
function of the temperature. 
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3 Results and discussion 



The gain against fiber lengths is calculated in the following way. Firstly, we take f(r) in 
Gaussian form, /(r) = exp{—r 2 /ojQ)/'wuj 2 where ujq is the spot size and the effective core 
area is tiuj^ = 33 jim 2 . Dopant distribution N(r) is also assumed to be Gaussian, N(r) ~ 
exp{—r 2 /uj 2 )/iiuj 2 . In addition, the ratio {uj/ujq) between Gaussian dopant distribution and 
transverse intensity profiles is selected as 0.3. Secondly, we obtain R and a s by using N(r) 
and f(r) for the relevant fiber parameters. Thus, the output pump power in Eq. (j2()j) is 
calculated with the different temperature values for the fiber length of 45 m. In this case, 
it is bear in mind that the ratio of cross-sections, which are belong to the signal beam, 
depends on the temperature. To calculate the parameter rj as a function of the temperature, 
we benefit by McCumber's theory, which gives a highly accurate relation between emission 
and absorption cross sections jllj . 

In the numerical calculations, we select the Al/P-silica erbium-doped fiber as an amplifier 
operated at the pump wavelength X p = 1480 ran and the input pump power P p (0) is fixed 
at 30 mW. The signal wavelength X s and the signal power P s (0) are taken as 1530 nm and 
10 ^W, respectively. The other parameters assigned to the fiber are given in Table 1 [T2*] . 
Moreover, we used the simulation programme OptiAmplifier 4-0 for generating P^ SE (L) only, 
and we set up the basic system seen in Figure 2 The energy difference between the 

sublevels of the metastable level (level 2) is assumed as 300 cm -1 in the room temperature 
for the simplicity. 

The results calculated for the various output pump powers and the intrinsic saturation 
powers as well as the parameters j3 and rj are given in Table 2. The variation of the signal 
gain against the fiber length is illustrated in Figure 3 for the temperatures —20 °C, 20 °C and 
60 °C. For a given pump and signal powers the gain decreases with increasing temperature 
within the range of L < 27 m. The difference between the maximum gains for —20 and 
60 °C is 0.67 dB. There is a temperature insensitivity for the length about L w 30 m for the 
relevant pump and signal powers. On the other hand, this temperature insensitive length is 
equivalent to the length at which the gain curves intersect each other (Figure 3). 
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4 Conclusion 



We have introduced a more accurate model including the temperature effect for the sig- 
nal gain of the erbium-doped fiber amplifier. In addition, we have shown the possibility 
of deriving an analytical solution of the propagation equations for some practical tempera- 
ture ranges. The temperature dependence of the output pump power is smaller than that 
of the intrinsic saturation power. Thus, in terms of the practical applications we can ne- 
glect the dependence of the output pump power on the temperature. However, it is taken 
into consideration that the gain performance of EDFAs strongly depends on the temperature. 
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Table 1: Typical fiber parameters for an Al/P-silica erbium-doped fiber (from Ref.|12|). 



Symbols Definations Values 





signal emission cross-section 


5.7xl(r 25 m 2 


< 


signal absorption cross-section 


6.6xl0~ 25 m 2 




pump emission cross-section 


0.87xl0" 25 m 2 


a" 
J p 


pump absorption cross-section 


2.44xl0" 25 m 2 


T 


life time 


10.8 ms 


N 


erbium concentration 


3.86xl0 24 m~ 3 




signal wavelength 


1530 nm 




pump wavelength 


1480 nm 




copropagating ASE power 


0.15 mW 
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Table 2: The relevant fiber parameters as a function of temperature. 



Temperature (°C) rj P P (L) P" 



- 20 0.306 0.845 2.308 mW 0.493 mW 

+ 20 0.357 0.862 2.311 mW 0.474 mW 

+ 60 0.406 0.879 2.314 mW 0.451 mW 
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Figure Captions 

Figure 1: Two level amplification system and main transition of erbium ion. 

Figure 2: Simulation setup for measurement of the co-propagating ASE power in an Er 3+ -doped 
optical fiber amplifier (from OptiAmplifier 4-0). 

Figure 3: Gain as a function of fiber length. P p (0) = 30 mW and P s (0) = 10 fj,W. 
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